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Abstract. We present Fabry-Perot and multi-object spectroscopy of the galaxies in Hickson compact group 31 (HCG31). 
Based upon our Ho- data cubes, galaxies A and C are a single entity, showing no discontinuity in their kinematics. Kinematically, 
galaxy E is probably a component of the A-rC complex; otherwise it is a recently detached fragment. Galaxy F appears, both 
kinematically and chemically, to have formed from material tidally removed from the A-l-C complex. Galaxies B and G are 
kinematically distinct from this complex. Galaxy Q also has a radial velocity compatible with group membership. Galaxies 
A, B, C, and F have nearly identical oxygen abundances, despite spanning a luminosity range of 5 mag. Galaxy B’s oxygen 
abundance is normal for its luminosity, while galaxy F’s abundance is that expected given its origin as a tidal fragment of the 
A-l-C complex. The oxygen abundances in galaxies A and C are also understandable if the A-t-C complex is a late-type spiral 
suffering strong gas inflow and star formation as a result of a tidal interaction. Given the kinematics of both the galaxies and 
the Hi gas, the oxygen abundances, and the position of galaxy G, we propose that an interaction of galaxy G with the A-l-C 
complex, rather than a merger of galaxies A and C, is a more complete explanation for the tidal features and other properties of 
HCG3L In this case, the A-tC complex need not be a merger in progress, though this is not ruled out. 

Key words, galaxies: abundances - galaxies: evolution - galaxies: formation - galaxies: individual: NGC 1741 - galaxies: 
clusters: individual: HCG31 


1. Introduction 


The Hickson compact group of galaxies 31 (HCG31; Hickson 


1982) has recently proven an interesting laboratory for study¬ 


ing various issues related to galaxy formation and merging. 
HCG31 is physically compact. Six of its members fall within 
the dimen sions of a typical gia nt galaxy, share a common H i 
envelope (Williams et al. 1991), and span a small range in ra¬ 
dial velocities ( |^ubin et al. 199C; |Tickson et al. 1992). It is 
therefore no surprise that HCG31 shows a wide array of ev¬ 
idence for galaxy interaction including: tidal tails and irregu¬ 
lar morphology, complex kinematics, and vigourous star for¬ 
mation, including starbursts with Wolf-Rayet spectral features 
(Rubin et al. 1990; Iglesias-Paramo & Vflchez [1997 ; Izotov & 


Thuan 1998; lohnson et al. 1999; lohnson & Conti 2000). 


Rubin et al. ( |1990| ) conclude that HCG31 will appear as a 
single normal galaxy within a few orbital periods. Here, we dif- 
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fer somewhat from this conclusion. Although we find abundant 
kinematical evidence for recent gravitational interaction among 
the group members, we suggest that there are other explana¬ 
tions besides imminent coalescence. Likewise, we find chem¬ 
ical abundances for the different components that are unlike 
those of normal dwarf irregulars, but again argue that their ori¬ 
gins may be easily understood. Given the abundant evidence 
for gravitational interaction, a merger may yet occur in HCG31, 
but we hnd that it need not be currently underway. 


In the following section, we present our data and their re¬ 
duction. In Sect. 3, we present our results. In Sect. 4, we discuss 
these results and their implications for the state of HCG31. In 
Sect. 5, we present our interpretation of the past and future evo¬ 
lution of HCG31. In Sect. 6, we summarize our conclusions. 


Throughout, we adopt the galaxy nomenclature of Rubin et al 


(199C) and we assume a distance of 58Mpc forHCG31, based 
upon a mean redshift of 0.0137 (Hickson et al. 1992) and a 
Hubble constant of 71 kms“* Mpc“‘ (Mould et al. 2000). A 
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Fig. 1. The field of HCG31 in the V-band. Galaxies A, B, C, 
D, E, and F are shown; galaxies G and Q are to the south¬ 
east and north, respectively, as indicated by the arrows. Deeper, 
higher-resolution images may be found in Rubin et al. (1990), 
Iglesias-Paramo & Vflchez (1997), Johnson et al. (1999), and 
Johnson & Conti (200C). The ellipses are meant to illustrate 
schematically the positions of galaxies A, C, E, F and the tidal 
dwarf candidates (denoted by “T”; Hunsberger et al. 1996). 
Objects within the component galaxies are labelled using the 
galaxy component letter and a sequence number, e.g., El and 
F2 are subcomponents of galaxy F. Where sub-components are 
labelled, the parent galaxies are not, to reduce the clutter. North 
is up and east is to the left, as indicated, and 1" = 281 pc. 
Unlabelled objects are foreground stars, except that north-east 
of galaxy D, which is a reflection. The filled circles indicate the 
locations of the apertures used for multi-object spectroscopy. 


preliminary version of some of our results was presented by 
Rosado et al. (1999). 


2. Observations and reductions 

The Fabry-Perot (FP) observations of HCG31 were carried out 
on 1 November 1997 at the f/7.5 Cassegrain focus of the 2.1 m 
telescope of the Observatorio Astronomico Nacional in San 
Pedro Martir, B.C., Mexico us ing the UNAM Scan ning Fabry- 
Perot Interferometer PUMA (Rosado et al. 1995). The scan¬ 


ning FP interferometer is based upon an ET-50 Queensgate 
Instruments etalon with a servo-stabilization system. The Ha 
line is observed in interference order 330, within a free spectral 
range of 19.89A, and sampled at 48 steps of 0.43A separation 
(18.9 km s '). The spectral resolution of these observations is 
38.4 km s“'. The effective finesse obtained with this setup was 
24. Since HCG31 subtends only 2' and we were careful to place 
it at the center of the instrument’s field of view (10') when ac¬ 
quiring the data cubes, the effective finesse of the observations 
is the same. In any case, PUMA does not show important varia¬ 
tions of the effective finesse across the field. A 30A interference 
filter, centered at the wavelength of redshifted Ha, was used to 
isolate the Ha emission line from HCG31. The detector was 


a 1024x1024, thinned Tektronix CCD. The image scale was 
0'.'59pixel *, yielding a 10' field of view. In order to increase 
the S/N of the observations, the detector was used with a 2 x 2 
pixel binning. 

We obtained two data cubes of HCG31 in Ha with ex¬ 
posure times of 48 minutes each (1 minute/channel). These 
data cubes were co-added to enhance the S/N of the faint re¬ 
gions. The seeing was about 1'.'2. The transparency conditions 
were rather good and we obtained both cubes during dark time. 
Nevertheless, we corrected both data cubes for transparency 
variations before co-adding them. This was done using two 
field stars (located outside the region shown in Fig. and we 
verified that the profiles were similar. 

We obtained wavelength calibration data cubes of a Ne 
lamp at the beginning and end of the observations, which also 
serve to check for possible equipment flexures. Since the red- 
shifted Ha emission of HCG31 differs from the wavelength 
of the calibration lamp, a phase shift correction was applied, 
which amounted to a shift in the zero-point of the velocities 


of 22 km s *. The CIGAFE software package ( Fe Coarer et al 


1993 ) was used to apply this phase shift correction, to remove 


cosmic rays, calibrate in wavelength, and construct the radial 
velocity cubes. We also used some routines from the Image 
Reduction and Analysis Facility (IRAF)||] for parts of the data 
reduction. 

The multi-object spectroscopy of the galaxies in HCG31 
was obtained with the 2.1m telescope of the Observatorio 
Astroflsico Guillermo Haro in Cananea, Mexico on 10 January 
1999. The FFOSC spectrograph was used, which is a trans¬ 
mission spectrograph employing a grism as the dispersing el¬ 
ement (Zickgraf et al. 1997). The spectrograph’s dispersion 
was approximately 5.5 A/pix and the spectra typically spanned 
from 4000A to 6600A, though the exact spectral range var¬ 
ied according to the object’s position within the spectrograph’s 
field of view. The detector was a 576 x 384 EEV CCD. 
Since the spectrograph has very poor response in the blue, 
no order-sorting filter was used. The objects were selected for 
spectroscopy using focal plane masks made from previously- 
acquired images. The only restriction for object selection is that 
they may not be aligned in declination since the dispersion axis 
is oriented east-west. Pairs of holes were cut for object and sky 
at the same right ascension so as to ensure identical spectral 
coverage. The holes were 3" in diameter. 

The objects selected for spe ctroscopy are show n in Fig. 
with the exception of galaxy Q ( |^ubin et~aL 1990 ), which is to 
the north of the field. The total integration time was 1.5 hours. 

The standard stars HD 19445, HD 74721, and HD 109995 


were observed for flux calibration. The standard stars were ob¬ 
served through masks cut for various object fields. 

The multi-object spectroscopy was reduced using the IRAF 
software package, specifically the noao.imred.specred pack¬ 
age. The bias images were first combined and the result sub¬ 
tracted from all of the images. Pixel-to-pixel variations were 


’ IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Universities for 
Research in Astronomy, Inc., under contract to the National Science 
Foundation. 
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Fig. 2. The velocity map of HCG31 in Ha. The velocities are 
heliocentric velocities. Panel (a) shows the velocities for galax¬ 
ies A, C, E, and F; panel (b) shows the velocities for galaxy B; 
and panel (c) shows the velocities for galaxy G. In all cases, 
the colour bars show the total range of radial velocities. The el¬ 
lipses are again meant to illustrate the positions of galaxies A, 
C, E, F, and the tidal candidates, as in Fig. [^. In general, both 
our radial velocities and velocity profiles are in good agree¬ 
ment with those of Rubin et al. (1990). The continuity of the 
kinematics for galaxies A, C, and E is particularly striking, and 
probably indicates that these three components are a single en¬ 
tity. Galaxy Q was not detected in our Fabry-Perot data cubes, 
while the Her line from galaxy D is shifte d outside our filter b y 
this galaxy’s much higher radial velocity ( Hickson et al. 1992 ). 
In all panels, north is up, east is to the left, and 1" =281 pc. 


then removed using spectra of the internal lamp. Subsequently, 
the sky was subtracted from each object. Spectra of the Ne- 
Ar lamp were used for wavelength calibration. Finally, the flux 
calibration was made using the standard star observations. 


3. Results 

Figs. ^ and Upresent the radial velocity and velocity dispersion 
maps resulting from our Fabry-Perot spectroscopy at Ha. To 
compute these maps, we computed the average intensity and its 
standard deviation over all channels at each spatial pixel. Then, 
we determined the median value of this standard deviation over 
all spatial pixels in the image and adopted this as representa¬ 
tive of the noise in the sky throughout the image. Finally, at 
each spatial pixel, the maximum of the intensities in all chan¬ 
nels was compared to this noise value. For those spatial pixels 
where the maximum intensity in some channel exceeded five 
times the noise value, a gaussian was fit to the velocity profile. 
The peak position and the standard deviation of this gaussian 
fit were adopted as the radial velocity and velocity dispersion, 
respectively, at that spatial position. 



Fig. 3. The velocity dispersion map of galaxies A and C in Ha. 
The colour scale shows the range in velocity dispersion values. 
Here, we characterize the velocity dispersion as the standard 
deviation of the gaussian fit to the velocity profile. Only galax¬ 
ies A and C showed significant internal variations in velocity 
dispersion, so they are the only ones shown here. The contours 
show the total Ha flux from the Fabry-Perot observations. As 
usual, north is up, east is to the left, and V — 281 pc. 


In Fig. m our radial velocities and velocity profiles are in 
good agreement with those found by ^^ubin et~^ (1990). In 
Fig. ^ we present several simulated profiles for the positions 
and position angles corresponding to Figs. 2 and 3 of Rubin 
et al. ( |1990| ). Galaxy Q was not detected in our data cubes, 
while galaxy D’s much higher radial velocity shifted its Ha 
line outside our filter. 

In Fig. H, we present the velocity dispersions in galaxies A 
and C. We present a map of the velocity dispersion for galaxies 
A and C only, since no other galaxies showed significant inter¬ 
nal variations in velocity dispersion. The velocity dispersions 
in galaxies B, E, and F are approximately 50kms \ similar 
to our instrumental resolution, while the velocity dispersion in 
galaxy G is slightly higher. 

Table lists the raw line intensities relative to Hfi and the 
raw radial velocities we measured in the H ii regions in the dif¬ 
ferent components of HCG31. The uncertainties quoted for the 
line intensities are Icr uncertainties and include the uncertainty 
in the fits to both the line and the reference line as well as the 
uncertainty in the continuum placement for both lines. The line 
intensities are not corrected for reddening. The names coincide 
with those indicated in Fig. Table also includes the redden¬ 
ing derived from the Balmer lines along with the abundance 
of doubly ionized oxygen and an estimate of the total oxygen 
abundance for those cases when we could derive an electron 
temperature. 

Before computing the reddenings listed in Table [^, we cor¬ 
rected the Balmer lines for underlying absorption according to 


,(d) = F„(T)x(l + i^)/(l + ^), 


( 1 ) 
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Table 1. Raw line intensities^* for the Hii regions in HCG31 measured from the multi-object spectroscopy 


A 


Al+Cl 

C2 

B2 

B3 

FI 

F2 

Q 

He 

3970 


12.51 ±0.80 






Hd 

4101 

10.44 ± 0.99 

18.42 ±0.86 



10.3 ± 2.7 



Hr 

4340 

29.2 ± 1.3 

38.4 ± 1.1 


33.6 ±3.3 

31.1 ±3.1 

31.9 ±5.3 


[Oiii] 

4363 

2.20 ± 0.93 

1.21 ±0.61 



5.1 ±2.5 



He I 

4472 

1.14 ±0.97 

2.70 ± 0.66 






HyS 

4861 

100.0 ± 2.8 

100.0 ± 2.9 

100 ± 16 

100.0 ± 3.5 

100.0 ± 2.6 

100.0 ± 4.5 

100 ± 32 

[Oiii] 

4959 

71.6 ±3.2 

90.0 ± 3.9 

135 ± 19 

117.2 ±4.6 

171.8 ±6.3 

142.9 ± 6.4 


[Oiii] 

5007 

218.1 ±5.4 

273.5 ± 6.9 

406 ± 47 

361.3 ±9.8 

521 ± 12 

429 ± 15 

157 ± 43 

He I 

5876 

15.68 ±0.43 

11.67 ±0.42 


12.0 ± 1.3 

13.83 ±0.91 

13.8 ± 1.6 


[Oi] 

6300 

8.77 ± 0.42 

4.77 ± 0.39 


10.7 ± 1.3 

8.06 ± 0.93 



[Oi] 

6364 

2.39 ± 0.28 

1.30 ±0.13 


1.51 ±0.92 

2.42 ± 0.67 



Ha 

6563 

473 ± 10 

326.4 ± 7.2 

625 ±71 


367.4 ±7.1 

385 ± 13 

470 ± 109 

[Nil] 

6583 

74.6 ± 3.9 

34.1 ±2.4 

57 ± 12 


11.8 ± 1.7 

12.9 ± 2.5 

51 ±20 

He I 

6678 

5.4 ± 1.0 

3.55 ± 0.34 



5.23 ± 0.88 



[Sii] 

6716 

47.8 ±2.1 

25.01 ±0.98 



24.8 ± 1.9 

23.6 ± 2.5 

146 ± 37 

[Sii] 

6731 

38.7 ± 1.8 

17.54 ± 0.77 



14.6 ± 1.7 

9.3 ±2.4 

125 ± 33 

VAHa) 

(km/s) 

3879 

3947 

4068 


3843 

3831 

4148 

E{B - V) 

(mag) 

0.42 ± 0.02 

0.10 ±0.02 

0.42 ±0.10 

0.68 ± 0.44 

0.21 ± 0.02 

0.26 ± 0.03 

0.11 ±0.22 

L([Olll]) 

(K) 

12500 ± 3200 

9100 ± 2000 



11900 ± 3200 



12 -t log(0++/H) 

7.56 ± 0.34 

8.12 ±0.39 



8.02 ± 0.38 



12 -t log(0/H)r. 

8.1 ±0.2 

8.3 ± 0.2 



8.1 ±0.2 




“ The line intensities are measured relative to HjS and are not corrected for interstellar reddening. 


where Fo(A) and Fc(A) are the observed and corrected fluxes 
and VTabs, Wa, and Wup are the equivalent widths of the under¬ 
lying stellar absorption, of the emission li ne in question, and 
of Hy6, respectively. We used Wabs = 1 -9A ( McCall et al. 1985 ) 
since our spectra had neither the resolution nor the signal-to- 
noise required to measure this directly. These corrected line 
intensities were then used to compute the reddening via 

log ^ = log 


m 

I(W) 


-0.4£(B-y)(Ai(d)-Ai(H/l)), 


( 2 ) 


where F(A) and /(A) are the observed and theoretical emission¬ 
line fluxes at wavelength d, E{B - V) is the reddening, and 
Ai(/l) = A(A)IE(B - V). Values of Ai(/l) were taken from the 
Schild ( 1977 ) reddening law. We used the Ho'/H/I ratio to de¬ 
termine the reddening, except when Ha fell outside our wave¬ 
length range, in which case the Hy/H/5 ratio was used. In all 
cases, we adopted theoretical ratios of - 2.86 and 

I(liy)/I(iy3) = 0.468, appropriate for electron t emperatures 
and densities of lO^K and 100cm“^, respectively (Osterbrock 


1989). 


The radial velocities given in Table are the raw radial ve¬ 
locities measured for the Ha line from the multi-object spec¬ 
tra, and are meant primarily to show that galaxy Q has a radial 
velocity consistent with group membership, and similar to the 
Hi velocity found by Williams et al. (1991). Galaxy Q is the 
faintest of the components we observed, which may explain 


why it is so faint in the Ha image presented by Rubin et al. 


(1990) and absent from our Ha data cubes. We cannot restrict 
galaxy Q’s morphological type since the V-band and Ha im¬ 
ages we obtained to identify H ii regions are of low spatial res¬ 
olution and signal-to-noise and were never calibrated. 


Table 2. Luminosities and oxygen abundances 


Galaxy 

B 

12 ± log(0/H) 

method** 

ref 

A 

15.62 ±0.5 

8.04 ± 0.06 

A 

1,3 

B 

15.35 ±0.2 

8.34 ± 0.20 

B 

1,4 

C 

13.44 ±0.5 

8.3 ± 0.2 

A,C 

1,5,6 

F 

18.2 ± 1.0“ 

8.1 ±0.2 

C 

2,6 


** The luminosity of component F is assumed equal to that of com¬ 
ponent E. 

’’ The oxygen abundances are obtained by the following methods: 
(A) measurements of [O ii] and [O iii] lines and the electron tem¬ 
perature; (B) measurements of [On] and [Oiii] lines only; (C) 
measurement of the [O in] lines and the electron temperature 


Reference s: (1) Hickson et al 1988 ; _(2)Jlglesi^- Paramo & Vflche 2 
1997t (3) |:zotov & Thuan||l99k (4) |Rubin et alj|l990|; (5) V acca~& 


Conti |l992t (6) Table 14 


The abundance of doubly ionized oxyge n is com puted ac¬ 
cording to the prescription given in |Tichei ( 1993 ) using the 
reddening-corrected line intensities and assuming an electron 
density of lOOcm^^. The total oxygen abund ances are esti¬ 
mated from Fig. 5 of ^obulnicky et al. ( 1999| ) in those cases 
when the electron temperature was measured. Generally, these 
total oxygen abundances are in good agreement with extant 
data (see Table ||). 

In Table ^ we list the luminosities and oxygen abundances 
for the component galaxies with measured oxygen abundances. 
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Fig. 4. This figure presents velocity profiles of galaxies A, B, 
and C extracted at the same positions and position angles as 
Rubin et al. (199C). Each panel indicates the position angle of 
the simulated slit. The agreement is very good for both the ve¬ 
locities and the velocity profiles, any discrepancies being of 
order lOkm/s (compare with their Figs. 2 and 3). We are not 
as sensitive to low surface brightness emission, so our profiles 
have a smaller spatial extent than those of Rubin et al. (1990). 


For each galaxy, we adopted what appeared to be the best 
abundance determination available, indicated in the footnotes 
to Table ^ using the precision quoted in the original papers as 
our guide. For galaxies A, B, and C, we adopted the Hickson 
et al. ( |1988[ ) B-band magnitudes at the 24.5 mag/□" isophote 
since these seem to be in better agreement with the photome- 



Mb (mag) 


Fig. 5. The relation between luminosity and metallicity for the 
galaxies in HCG31. The components are labelled. The solid 
line is the relation for dwarf irregulars found by Richer & 
McCall (1995), while the dashed line is an extrapolation of that 
relation. The galaxies in HCG31 are clearly discrepant, with 
components A and C much too luminous for their abundances. 


try of Iglesias-Paramo & Vilchez (1997). No photometry ex¬ 
ists for galaxy F. Here, we assume it has the same luminosity 
as galaxy E (Iglesias-Paramo & Vilchez 1997), which, though 
it is a rough approximation, is sufficient for our purposes. In 
Table ^ we adopted the uncertainties in the luminosities of 
galaxies B and C from Hickson et al. (1988). For galaxy A, 
we set the uncertainty in luminosity equal to that for galaxy 
C due to the difficulty of disentangling the two galaxies. We 
arbitrarily estimate an uncertainty of 1.0 mag in the luminos¬ 
ity of galaxy F. In Fig. ^ we plot the oxygen abundance as a 
function of the galaxy luminosity. We converted from apparent 
to absolute magnitudes using an apparent distance modulus of 
34.02 mag, which accounts for our adopted distance and an ex- 
t inction of = 0.21 mag, b ased upon E{B - V ) - 0.055 mag 
( Schlegel et al. 1998 ) and the Fitzpatrick ( 199S ) reddening law 
{Rv = 3.1 and E{B - V) = 1 mag). We also show the relation 
between oxygen abundance and luminos ity for dwarf irregulars 
for comparison ( Richer & McCalj 1995 ). 


4. The state of HCG31 today 


Our most surprising result is that galaxies A and C appear to be 
a single kinematical entity. In Fig. ^ their kinematics show no 
discontinuity. We interpret this as evidence that galaxies A and 
C are a sin gle entity at present. Mendes de Oliveira & Amram 
( 200C ) and Amram & Mendes de Oliveira ( 2000 ) find similar 


results from similar kinematical data. 

Another surprise is that galaxy E’s position and velocity 
are compatible with it also being an integral part of the A-tC 
complex. If galaxy E is in fact a separate entity, it is just now 
separating from the A-tC complex. In any case, its kinematics 
confirm that the extension of the A-tC complex towards galax¬ 
ies E and F is a tidal tail. 
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Like galaxy E, the tidal dwarfs to the north-east of galax¬ 
ies A and C identified by Hunsberger et al. (1996) also have 
kinematics consistent with their being integral components of 
the A-hC complex. Furthermore, given the red continuum im¬ 
age presented by Johnson & Conti (2000), it is not evident that 
at least the nearer of these supposed tidal dwarfs are separate 
entities. 

Meanwhile, galaxy F is an excellent example of a tidal frag¬ 
ment. Considering its position and radial velocity (Figs. [| and 
it seems clear that galaxy F is a tidal fragment of the A-hC 
complex that detached some time ago. Remarkably, galaxy F 
has the same oxygen abundance as galaxies A and C, despite 
being more than 5 mag fainter than galaxy C (cf. Table || and 
Fig. ||). In addition, since an underlying stellar population has 
yet to be detected in galaxy F (Johnson & Conti 2000| ), it can¬ 
not be responsible for its own chemical enrichment. All of these 
observations argue that galaxy F is a tidal fragment of the A-hC 
complex. 

Galaxy B exhibits solid-body rotation. Given its broadband 
morphology (j^ubin et al. 1990; Iglesias-Paramo & Vflchez 


Galaxy B, meanwhile, appears more or less normal. It has 
a B-band luminosity at the bright end of the range found for 
dwarf irregulars, about 0.8 mag brighter than the LMC’s, and 
an oxygen abundance similar to that of the LMC. Its Ha lu¬ 
minosity, 1.0 X 10"^' ergs^' (Iglesias-Paramo & Vflchez 1997, 
corrected to our adopted distance), implies a global star for¬ 
mation rate of 0.8 Mq yr"' (Kennicutt||1998|), which is not un¬ 


usual for a late-type galaxy (Kennicutt 1983). However, its 
colours are significantly bluer than those of typical dwarf ir¬ 
regulars or late-type spirals (Iglesias-Paramo & Vflchez 1997 


|de Vaucouleurs et al. 1981; Roberts & Haynes 1994), being 
a near-perfect match to those of a single stellar population of 
metallicity Z = 0.008 and an age of 40Myr (Bressan et al 


1994). Thus, galaxy B is currently forming stars more actively 


than it has formed them on average in the past, perhaps as a 
result of the group environment. In spite of this, galaxy B’s re¬ 
cent history of star formation appears to be distinct fro m that of 
the A-i-C complex ( Iglesias-Paramo & Vflchez 1997 ; Johnson 
& Conti 2000). 


1997), it is likely a dwarf spiral or irregular seen nearly edge- 5. The past and future 


on. Kinematically, galaxy B is distinct from the A-i-C complex 
and is counter-rotating with respect to it. 

Galaxy G appears to be a small disk galaxy seen nearly 
face-on since the Ha velocity has only a very small gradient 
across its disk. 

The velocity dispersion map (Fig. |]) of the ionized gas is 
more difficult to interpret. The highest velocity dispersions oc¬ 
cur along a ridge between the main bodies of galaxies A and C. 
Whether this higher velocity dispersion is due to a more com¬ 
plicated gas distribution or kinematics is indeterminate from 
our observations. 

It is evident from Fig. ^ that the oxygen abundances in 
the galaxies in HCG31 are unusual. Although they span 5 mag 
in luminosity, all of the galaxies have nearly identical oxygen 
abundances. This conclusion is independent of the uncertainty 
in separating the luminosities of galaxies A and C, since, were 
we to consider their combined luminosity and either of the 
oxygen abundances observed (Table ||), the resulting position 
would remain highly discrepant compared to the relation for 
normal dwarf irregular galaxies. However, galaxy C’s luminos¬ 
ity is much greater than that of most dwarfs and more like that 
of a typical late-type spiral. If the A-hC complex were a late- 
type spiral, its oxygen abundance may be more easily under¬ 
stood. Given the strong inflows and consequent high star for¬ 
mation rate expected due to tidal interactions (e.g., Mihos & 
Hernquist 1996| ; Barnes & Hernquist 1996), the mass of gas 
expected in the outer disks of late-type spirals, and the low 


The kinematical evidence argues strongly that the A-i-C com¬ 
plex has undergone very significant, recent, tidal interaction, 
but that it is now a single entity. In addition, the kinematical 
properties of galaxies E, F, and the tidal candidates identified 
by Hunsberger et al. (1996) all provide evidence of their origin 
within the A-i-C complex. What, then, are the possible origins 
of these tidal interactions, and might these shed some light on 
the nature of HCG31 ? A related issue is that of the nature of 
the original components of HCG31. 

The global H i mass and blue luminosity of HCG31 provide 
some clues regarding the original composition of HCG31. The 
Hi mass of HCG31 is s Rikingly high, 2.5 x IO'^Mq adjusted 
to our adopted distance ( Williams et al. 1991), which would b e 
typical of a massive Sb or Sc spiral ( Roberts & Haynes 1994 ). 
The oxygen abundances, however, are incompatible with those 
of such mas sive galaxies save those found in their ext reme 
outer disks ( Garnett et al. 1997 ; l^erguson et al. 1998| ; van 
Zee et al. 1998 ). The blue luminosities of galaxies A-i-C, B, 
and G im ply m orphological types Scd, Sm, and Sd (Roberts 
& Haynes 1994 ), though these are all undoubtedly affected by 
their significant recent star formation. The total H i mass is also 
consistent with the sum of three late-type spirals (or one such 
spiral and two bright irregulars; Roberts & Haynes 1994 ). As 
noted above, the oxygen abundances of Ah-C and B are com¬ 
patible with such progenitors, though that for A-hC would be 


oxygen abundances observed in these regions (Ferguson et al. 


at the low end of t he range expected ( Richer & McCall 1995 
[Garnett et~d! 1997). Finally, the global LbIM{}\ i) for the entire 


1998; van Zee et al. 1998), the unusual positions of galaxies 
A and C in Fig. ^ are understandable as the result of dilution 
due to strong gas inflow and the transient high luminosity due 
to a starburst. At any rate, galaxy C’s position is compatible 
with the extreme low end of the abundance range for late-type 


spirals (Garnett et al. 1997). Given galaxy F’s origin within the 


A-I-C complex, it should naturally have the same oxygen abun¬ 
dance as the A-I-C complex. 


group is approximately unity, a value typical of dwarf irregu¬ 
lars (Roberts & Haynes 1994). Despite the uncertainty due to 
the recent star formation, all of the available evidence points 
toward the group being initially composed of late-type spirals 
and irregulars. 

The kinematical evolution of HCG31 is a more ambiguous 
puzzle. Based upon the galaxy morphologies and less complete 
kinematical data than we present here, an on-going merger be¬ 
tween galaxies A and C has been invoked to explain the tidal 
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features observed (e.g., |^ubin et"^ 199C). However, this ex¬ 
planation leaves a number of issues unresolved, the two most 
important of which are the coincidence of galaxy G with a 
peak in the Hi distribution at the end of a tidal tail and the 
very similar recent star formation histories in galaxies A, C, 


ity dispersions we measure in these galaxies (Fig. ^ Sect. 3), 
and it is not obvious whether this encounter velocity should 


provoke such a strong interaction (e.g., Binney & Tremaine 
|1987). The time scale could be lengthened, and the transverse 


E, F, and G (Williams et al. 1991; Iglesias-Paramo & Vflchez 


velocity lowered, if there is a delay between closest approach 
and the initiation of star formation, as some models indicate 


1997; Johnson et al. 1999). The continuity of the kinematics (e.g., Mihos et al. 1992; Mihos & Hernquist 1996). On the 


of the A-hC complex is an additional important complication 
since it implies that the interaction would be in an advanced, 
settled stage, in contrast to numerical simulations that indicate 
that most of the star formation ought to occur at the onset of 
the interaction (Barnes & Hernquist 1996; [Vlihos & Hernqui^ 


1996). Likewise, the blue colours of the star clusters in galax- 


i es A and C imply th at any merger should be in an initial phase 
( Johnson et al. [l9^ ). 

In vi ew of the continuous kinematics of t he A-i-C complex 
(see also 4mram & Mendes de Oliveira 2000 ), we propose that 
an encounter between it and galaxy G is a more likely origin for 
the tidal interaction that is so evident. This encounter scenario 
naturally accounts for the position of galaxy G, embedded in 
the south-east tidal tail of the H i distribution ( Williams et al. 


spite their physical separation of up to 40kpc, since such an 
encounter would provoke strong gas inflow and presumably 
vigourous star formation in galaxy G, the A-hC complex, and in 
the tidal material (Mihos et al. 1992; Mihos & Hernquist 1996; 
Barnes & Hernquist 1996). Such strong gas inflow would also 
explain the low oxygen abundance measured in the A-hC com¬ 
plex as a result of dilution by metal-poor gas. The coordinated 
star formation, the global morphology and kinematics of the 


Hi gas (Williams et al. 1991), and the positions of galaxies 
Ah-C, E, F, and G are all very reminiscent of the fly-by and 


( |1996| ) , and [Barnes & Hernquist| ( |1996| ). Finally, O’Halloran 


merger simulations of Mihos et al. ( 1992 ), |Mihos & Hernquist 


other hand, neither the velocity nor the spatial extension be¬ 
tween galaxy G and the A-i-C complex are unusual. HCG26 
is very similar (Williams & van Gorkom 1995), while HCG16 
has a similar morphology on a larger scale (Verdes-Montenegro 


et al. pOOip . NGC 7318b in Stephan’s Quintet (HCG92) has 
been stripped of its Hi despite a velocity of order 1000kms~* 


relative to the other group members (Sulentic et al. ZOOl). 


Likewise, UGC 7636 in the Virgo cluster has had its interstel¬ 
lar medium stripped by NGC 4472 though their relative ve¬ 
locities differ by about 700kms~* (Patterson & Thuan 1992 


jMcNamara et al. 1994 ; Lee et al. 200C). Therefore, although 


1991). An encounter between galaxy G and the A-hC complex 


also easily explains the very similar histories of the most recent 
star formation in galaxies A-hC, E, F, and G noted by Iglesias- 
Paramo & Vflchez (1997) and lohnson & Conti (200C), de- 


the response of the H i in fly-by encounters appears to be rela¬ 
tively unexplored theoretically, observational examples appear 
reasonably common, lending credence to the idea that an inter¬ 
action between galaxies G and Ah-C could be the origin of the 
tidal debris seen in HCG31. 

If the gravitational interaction between galaxy G and the 
Ah-C complex is the origin of the tidal debris in HCG31, it is 
then also possible that the A-hC complex has always been a 
single entity and need not represent the fusion of two galaxies, 


partially confirming the hypothesis of Plana et al. (200C). The 


simulations of fly-by and merger interactions indicate that very 
strong gas inflows may be induced (|Mihos et ^ 1992; Mihos 


are indeed observed in HCG31 (Williams et al. 1991), though 


& Hernquist |1996tparnes & Hernquist 1996) and such inflows 


et al. (2002) conclude that the burst of star formation in the 


Ah-C complex may be evolving to a post-starburst phase, based 
upon ISO data, a result compatible with the encounter sce¬ 
nario if galaxy G is now moving away from the Ah-C com¬ 
plex. Therefore, unless galaxy G occupies its position purely 
by chance, we propose that an encounter between it and the 
Ah-C complex accounts for more of the properties of HCG31 
in a more natural fashion than does a merger of galaxies A and 
C. 

An important uncertainty in the encounter scenario is the 
transverse velocity implied between galaxy G and the Ah-C 
complex. Given their very similar radial velocities, the bulk 
of their relative motion must be in the plane of the sky (Fig. 


the HI kinematics are slightly compromised since their spatial 
resolution is insufficient to clearly separate the gas belonging 
to galaxy B from that belonging to the Ah-C complex. The dis¬ 
tribution of the star formation observed in the Ah-C complex 
could simply be a typical result of the gas inflows predicted 
by the numerical simulations. Likewise, the low oxygen abun¬ 
dances could easily result from dilution due to the inflow of 
metal-poor material from the outer disk. 

What might be the eventual outcome of HCG317 If an en¬ 
counter between galaxy G and the Ah-C complex is the origin of 
the tidal interactions, and if their orbits are unbound, the Ah-C 
complex will likely evolve similar to the simulations of ring 


galaxies (e.g., Struck 1997), and the Ah-C complex should re¬ 


main as a disk galaxy. If, on the contrary, their orbits are bound, 
a merger that will likely also consume galaxies E, F, and the 


tidal dwarfs postulated by Hunsberger et al. (1996) is the prob¬ 


able result. Finally, if the interaction between galaxies A and C 
is the origin of the tidal debris in HCG31, as has been usually 


Rubin et al. 

1990; ^ickson et al. 

1992 

). At a distance supposed (e.g., 

Rubin et al. 

1990; 

lohnson et al. 

1999 


of 58Mpc, the 2^4 separation of galaxy G and the Ah-C 
complex corresponds to a linear distance of 40kpc. If the 
ages of the most recent star formation, < 10* years (Iglesias- 


Paramo & Vflchez 1997; Johnson et al. 1999), are taken as the 


time elapsed since the putative encounter, the separation im¬ 
plies a transverse velocity in the plane of the sky exceeding 
400kms“*. This velocity greatly exceeds the internal veloc- 


well underway. 

In either case of a merger, we expect that a low mass ellip¬ 
tical would be the eventual outcome, as suggested by Rubin 
et al. (1990|). The group’s Hg luminosity is 1 x 10^^ergs 


(Iglesias-Paramo & Vflchez 1997, adjusted to our distance), 
implying a star formation rate of 7.9 Mq yr“* (ECennicutt 1998). 
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Although high, this rate of star formation is sustainable for a 
long time, at least another 3 Gyr given the group’s H i mass of 
2.5 X 1O'°M0 (Williams et al, 1991). Most of this Hi mass 


is contained within an area of 3' x 1' (Williams et al, 1991), 
8.6 X 10^ pc^ at our adopted distance, implying a gas surface 
mass density of about 30 Mq pc“^. Supposing a conversion ef- 
hciency from Hi gas into stars of approximately 10% at this 
surface density (ECennicutt 1998) and that at least 50-60% of the 
mass converted to stars will produce long-lived stars (Koppen 
& Arimoto 1991), we infer that the eventual outcome could 
be a galaxy with a stellar mass of several 10® Mq. This mass is 
several times larger than the current mass estimates: the colours 
from Iglesias-Paramo & Vflchez (1997) are completely domi¬ 
nated by the light from young stars, and therefore imply only 
lower limits to the stellar mass of order 3x10*^ Mq, based upon 


the colours of the single stellar populations from Bressan et al. 


( 1994 ). Therefore, it is reasonable to expect an eventual metal 
production at least equal to that which has already occurred, 
accompanied by a similar increase in the chemical abundances. 
Several Gyr after the cessation of star formation, a galaxy with 
a stellar mass of sev eral 10® Mq would have a luminosity of 
Mb ~ -18 mag (e.g., Bressan et al. 1994) and its oxygen abun¬ 
dance would place it above the metallicity-luminosity relation 
for dwarf irregulars. Given that the velocity gradients in the H i 
gas found by Williams et al. (1991) imply generally radial in¬ 
flow, rather than rotation, we infer that the probable outcome 
will be an elliptical galaxy. Although the foregoing is entirely 
speculative, it illustrates that HCG31 could very plausibly pro¬ 
duce a typical low-luminosity elliptical with normal chemical 
abundances (e.g., |Ticher et'd!||1998[ ). 


6. Conclusions 


Somewhat at variance with previous studies, it is less evident 
to us that HCG31 is a clear example of an on-going merger. 
It is undeniable, however, that there is signihcant tidal inter¬ 
action underway among the members of this compact group. 
Given their kinematics, gal axies A, C, and the tida l dwarfs to 
the north-east identihed by Hunsberger et al. ( 1996 ) all appear 
to be a single entity at present (the A-i-C complex). Galaxy E 
is probably also an integral part of this complex, otherwise it is 
a tidal fragment that has very recently separated. Galaxy F has 
both kinematics and chemical abundances that indicate that it 
is a tidal fragment that was part of the A-i-C complex in the 
relatively recent past. At this point, galaxies B and G remain 
kinematically distinct fro m the A-iC compl ex. Finally, we find 
a velocity for galaxy Q ( |^ubin et~n! 199C| ) that is compatible 
with it also being a member of HCG31. 

We propose that a gravitational interaction between galaxy 
G and the A-hC complex is the most complete explanation of 
the events occurring in HCG31. This scenario naturally ac¬ 
counts for the positions and kinematics of galaxies A, C, E, F, 
G, and the tidal dwarf candidates of Hunsberger et al. (1996), 


the morphology and kinematics of the Hi (Williams et al. 


1991), the chemical abundances, and the coordinated star for¬ 


mation in galaxies separated by up to 40 kpc (Iglesias-Paramo 
& Vflchez 1997; lohnson & Contj 200C). In this case, the A-hC 
complex may have always been a single entity, a late-type spi¬ 


ral, and its current irregular morphology, high luminosity, blue 
colours, and anomalous chemical abundance may simply be a 
consequence of multiple bursts of star formation as a result of 
strong gas inflows induced by its interaction with galaxy G 
(Barnes & Hernquist 1996| ; |VIihos & Hernquist 1996| ). If this 
is the case, a merger may still eventually occur between galaxy 
G and the A-i-C complex given the ample evidence for a strong 
gravitational interaction. While there is no evidence that the 
kinematical properties of HCG31 are not the result of an ongo¬ 
ing merger between galaxies A and C, many additional details 
are not accounted for by this solution, particularly the coordi¬ 
nated star formation and the location of galaxy G. 

Based upon its oxygen abundance, its broad-band and Her 
luminosities, and its kinematics, at present galaxy B appears 
to be a relatively normal late-type spiral or bright dwarf irreg¬ 
ular seen nearly edge-on. Its evolution is probably now being 
affected by the group environment, since its colours are very 
blue and indicate that its current rate of star formation exceeds 
its historical average rate. 
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